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Abstract 

We evaluate the p-wave K~N amplitudes from the chiral Lagrangians and from 
there construct the p-wave part of the K~ nucleus optical potential plus a small 
s-wave part induced from the elementary p-wave amplitude and the nuclear Fermi 
motion. Simultaneously, the momentum and energy dependence of the s-wave optical 
potential, previously developed, are taken into account and shown to generate a small 
p-wave correction to the optical potential. All the corrections considered are small 
compared to the leading s-wave potential, and lead to changes in the shifts and 
widths which are smaller than the experimental errors. A thorough study of the 
threshold region and low densities is conducted, revealing mathematical problems 
for which a physical solution is given. 

PACS: 13.75.Jz; 21.30.Fe; 21.65.+1; 25.80.Nv; 36.10. Gv 
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1 Introduction 



The problem of kaonic atoms has attracted considerable attention along the years fl|-|9| and 
has regained interest recently due to the new perspective that the use of chiral Lagrangians 
has brought to the problem of the kaon interaction with a nuclear medium |10|-[T3|. The 
need to obtain accurate results for the kaon selfenergy in a nuclear medium in view of 
the possibility to get kaon condensates in neutron-proton stars 



renewed interest in the subject. Similarly, the interpretation of the enhancement of 

the K~ yields in heavy ion reactions [19,^0| relies on the value of the K~ selfenergy in the 
nuclear medium. One of the characteristics of the KN interaction at low energies is the 
strong dominance of the s-wave amplitude, and the KN cross section to different channels 
can be very well described with only s-waves up to momenta of the order of 200 MeV or 
more [pl|-|2~5* . 

The dominance of the s-wave in the elementary interaction has been the justification 
for using traditionally s-wave K nucleus optical potentials fl]-|9[], by means of which good 
agreement with data can be obtained. This is in contrast with the situation in pionic atoms, 
since the relative small contribution of the s-wave part of the potential, together with a 
sizeable p-wave part, make this latter contribution quite important in the interpretation 



of the pionic atom data |f26|-|28| . 

Theoretical evaluations of the optical potential for a particle in a nucleus usually start 
from the impulse approximation, or tp form of the potential, with t the elementary scat- 
tering amplitude of the particle with the nucleons averaged over isospin and the Fermi 
motion of the nucleons. Yet, it is known that the explicit consideration of the Fermi mo- 
tion leads to corrections to this result. Certainly there are other higher order corrections in 
the density from mechanisms involving many nucleons, which one can systematically tackle 
using many body techniques |28] , but even at the level of one single scattering the explicit 
consideration of Fermi motion brings corrections to the optical potential beyond the tp 
approach. A detailed study of these correction terms for the case of pionic atoms was done 
in P5| , where it was found that the s-wave elementary amplitude induces a p-wave term 
in the optical potential and, similarly, the p-wave elementary amplitude induces an s-wave 
term. Yet, these corrections are small compared to the original impulse approximation and 
even small compared to the p 2 terms which account for pion absorption. 

The aims of the present paper are to investigate the nonlocal effects induced form the 
s-wave part of the K~ nucleus optical potential and also to evaluate the p-wave part 
of it. Firstly, we begin from the part of the potential evaluated from the s-wave K~N 
interaction in and derive from there the nonlocal corrections originated from the u) 
and k dependence of this potential. These non-local terms were evaluated in |Hj and there 
they were found to be large and quantitatively important in the analysis of the spectra of 
kaonic atoms. We also evaluate the part of the potential originating from the elementary 
p-wave K~N interaction and deduce from there the p-wave term of the optical potential 
plus an induced s-wave part due to the Fermi motion of the nucleons. 

We shall see that all these corrections are rather small and do not appreciably change 
the results obtained from the s-wave part of the potential alone. 
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One important finding is that, in order to properly evaluate these corrections, one must 
pay a special attention to the region of low densities. Consequently we discuss in detail the 
problems that one faces implementing the low density limit around a threshold, particularly 
in the derivatives of the optical potential. 

Finally, we study the effects of the inclusion of the £*— h excitation in the self-consistent 
calculation of the s-wave K~ potential. We find that, though best-fits of better quality 
can be obtained, this new theoretical potential leads to an overall acceptable description 
of the measured shifts and widths. 



Brief summary of the KN and K nucleus interac- 
tions 



For the elementary KN interaction we follow the chiral unitary model of |24 ]. This work 



follows closely the steps of [23|, where the usefulness of combining unitarity in coupled 



channels with the chiral Lagrangian was made manifest, but uses an enlarged basis of 
coupled channels. The inclusion of all the channels formed by the octets of the pseudoscalar 
mesons and stable baryons made it possible to obtain good solutions in |3] by means of only 
the lowest order chiral Lagrangian and a suitable cut off to regularize the loop integrals. On 
the other hand, in |23[ the effect of the unopened channels was accounted for by including 



higher order terms in the chiral Lagrangian. 

The K~N t matrix, ty, is obtained in [|23J] through the iteration of the Lagrangian in 
a coupled channel Bethe-Salpeter equation. 

The K~ selfenergy is evaluated in detail in ref. |13[ for nuclear matter by means of the 
integral 

n(M = 25:/ -j^p&p) t{ K- N (k,p,pl (i) 

where N stands for protons or neutrons and t^ m ' is the K~N scattering matrix in the 
nuclear medium with density p. In Eq. ([!]), n(p, p) denotes the occupation probability 
of momentum states in the Fermi sea in the nuclear medium at finite density p. Only 
the s-wave amplitude is considered in this evaluation. In section 5 we shall work out the 
contribution from the p-wave interaction. 

The t( m ) matrix is evaluated from the Bethe-Salpeter equation, but modifications are 
done in the meson and baryon propagators of the loops to incorporate the medium effects. 
The states allowed in the loops are KN, 7rS, ttA, 77E, 77A, KE,. 

The medium modifications implemented are the following: 
1) Pauli blocking in the nucleon propagators. This effect was proved to be very important 



in |10| ) [TT| . Indeed, forcing the intermediate nucleon states to be on top of a Fermi sea costs 
more energy, and the net effect is a shift to higher energies of both the real and imaginary 
parts of the K~p amplitude which is dominated by the A(1405) resonance below threshold. 
The shift of the real part automatically produces an attractive K~ selfenergy already at 
very small densities. 
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2) However, if the K~ selfenergy acquires a negative value then it costs less energy to 
produce the A (1405) resonance hence producing a shift of the K~p amplitude toward 
lower energies. A selfconsistent evaluation becomes then necessary as shown in [12], where 
it was found that the consideration of the K~ selfenergy together with Pauli blocking on 
the nucleons left the position of the A (1405) resonance basically unchanged. In the 
K~ selfenergy is also considered and a selfconsistent evaluation is also done. 

3) In addition to the former ingredients new effects are considered in [i~3|, i.e. the pion 
selfenergy in the 7rE, 7rA channels is also taken into account allowing the pions to excite ph, 
Ah and 2p2h components. Furthermore, the difference of binding between the nucleons or 
E and A hyperons is also incorporated. 

The results obtained are qualitatively similar to those found in |12| except that the 
imaginary part of the K~p amplitude becomes even wider and essentially flattens at full 
nuclear density p ~ p . 



3 Nonlocal terms associated to the s— wave part of the 
potential 

As described in the former section the optical potential for the K~ nucleus interaction was 



evaluated in nuclear matter as a function of the density in |13| and, in order to apply it to 
finite nuclei, the local density approach was used in ||, something justified for the s-wave 
potential as discussed in [P^j . The K~ selfenergy obtained in (JTSJ from the interaction of 
the K~ with protons and neutrons in symmetric nuclear matter has an explicit dependence 
on k° = uj and k, the energy and momentum of the antikaon. However, in order to solve 
the Klein-Gordon equation (KGE) to obtain energies and widths of the kaonic atoms in || , 
the potential was evaluated at the K~ threshold (u = rriK,k = 0). In what follows we 
derive the corrections to the optical potential from the consideration of the explicit u and 
k dependence of the kaon selfenergy in the nuclear medium. 
We write the K~ selfenergy in nuclear matter, n, as 



U(u, k, p) = 2ujV opt = Tl{m K , 0, p) + b(p) k 2 + c(p) (u - m K ) 



with 



Kp) 



dk 



c{p) 



(w=mj{, k=0) 



dU 



(2) 



(3) 



(u>=mj(, k=0) 



where the second order corrections in k 2 and (u — tuk) have been neglected. 

— # 

Once at this point the momentum k is not defined for the bound K~ in the atom and 
instead it becomes an operator. In detailed studies of finite nuclei one can trace the origin 
of this operator and how it acts on the density-dependent functions of the potential or the 
kaonic wave function. For instance, in |28| one can see that the k 2 which appears in the 
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p-wave part of the 7r-nucleus optical potential evaluated in nuclear matter corresponds in 
the equivalent finite nucleus calculation to the combination 



k 2 CM f(p) 



-Vf{p{r))V+ l -eV 2 fW)) 



m-n 



M, 



(4) 



N 



In the present case the evaluation of the K~ selfenergy in finite nuclei, with all the 
effects considered in [|l"2"|, |13l| plus the requirement of selfconsistency is a rather involved 
task, which would become advisable should these nonlocal effects be too big. Yet, as we 
shall see, the effects are small, smaller than present experimental uncertainties in the data, 
and thus the estimates which we shall perform here are sufficient to establish the relevance 
of these effects. One of the handicaps of having evaluated the selfenergy in infinite nuclear 
matter is that we do not know to which kind of VV operator will the factor k 2 correspond 
to. In order to estimate the size of the corrections and the uncertainties, we shall work 
with some different assumptions which are guided by the results obtained for pionic atoms 
in the translation from infinite matter to finite nuclei. The types of operators used for the 
b k 2 term of the .^"-selfenergy in Eq. (0) are shown below 



(5) 



a)bk 2 - 


* -V&V 


b) b k 2 - 


+ -V&V- i(A6) 


c) b k 2 - 


+ -bV 2 



The form a) is the Kisslinger type of interaction k ■ bk, the form b) appears in [2{J by 
using the Wigner transform of the symmetrized form (bk 2 + k 2 b)/2, the form c) allows the 
k 2 operator to act directly on the kaonic wave function, and thus has a special physical 
significance. 

We can also adopt a different point of view and, since the results of Eq. flip already 
come from using local approximations implicit in the use of a local Fermi sea, we can 
convert the nonlocal potential of k 2 into a local, but energy dependent, potential. The 
trade between nonlocal and local energy dependent potentials is a technique often used in 
many body theory It is based on the use of the Schrodinger equation (KGE in our 
case). The KGE is written here as 



[-V 2 + fi 2 + n(r)]0(f ) = [uj- V c (r)] 2 <f)(r ) . (6) 

where \i is the kaon-nucleus reduced mass and Vc(r) is the Coulomb potential with a finite 
nuclear size and vacuum-polarization corrections. 

In view of Eq. (^J) , and the fact that the meaning of k 2 in Eq. (|2]) is the kaon momentum 
squared, we can take for it the expectation value of —V 2 for the kaonic wave functions. 
This leads to 



d) b k' 



n 



(7) 
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Also since the potential is complex we can take the more symmetrical situation (bk 2 + 
k 2 b)/2, in which case we obtain the real part of the expression in Eq. (|7|) 



e) b k' 



b Re 



The energies and shifts for the energy dependent potential can be obtained by iteration. 
Given the smallness of the pieces under consideration the convergence is extremely fast. 



4 Threshold behaviour and the low density limit 

Eqs. (|3|) require the evaluation of the derivatives of the kaon selfenergy with respect to k 2 
and k° at threshold. However, these magnitudes are problematic when one goes to small 
densities, as we shall see. The problems stem from the behaviour of the free t matrix 
at threshold which has a cusp. This is related to the contribution of the elastic channel 
K~N — * K~N and hence, for the discussion, one can neglect all the other channels. 
Also the discussion here is completely general and thus we do not particularize to the 
K~N — > K~N reaction. 



4.1 Behaviour of the elastic rescattering terms 

To simplify the discussions we shall take unity for the transition potential K~ N — > K~ N . 
We begin with the study of the second order rescattering term with this potential, see 
Fig. [|. This is the term studied in the appendix of Ref . [30 , which there led to coefficients 



b and c behaving like k\ in the limit of low densities, a result that would seem to violate 
the low density theorem. 

In order to clarify the problems let us write the kaon selfenergy for this rescattering 
term 

U res (k° k p)=4 f ^ I 1 n(p,p)[l-n(p + k-q,p)] 

' J (2ixfJ (2tt) 3 2u(v) k° - v(q) + e(p) -e{p + k-q) +ir) ' 

where we use relativistic energies for the meson but non-relativistic ones for the nucleon, 
in order to be able to use explicit formulae for the Lindhard functions. In the following, 
we will denote the meson and nucleon masses by m and M, respectively. We have written 
explicitly the integral over the Fermi sea momentum, p, and over the meson momentum, 
q. Let us examine the behaviour of the q integral in Eq. @ in the limit of zero Fermi 
momentum and k = 0. In such a case one gets 

t rea (k°, k = 0) = ——r / q 2 dq—— 5 , (10) 

y 2^(g)fcO-o;(g)-^ + zr / 1 ' 

which would give us t res , which is the contribution of the rescattering term to the free t 
matrix, as a function of k°. Now let us find out the derivative of this function with respect 
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Figure 1: Second order rescattering contribution to the meson selfenergy. 



to k° and particularize for k° = m. Commuting the partial derivative and the integral, 
and concentrating in the region of small q, one finds an infrared divergence of the type 



Q 

Certainly, as one puts a finite, even if small, Fermi momentum, the divergence disap- 
pears, and one finds a 1/kp behaviour, similarly to the derivative of the function I(u, 0) 



with respect to u in the appendix of |30|. Obviously, the divergence is recovered in the 
strict limit of kp = 0. This is a subtle problem. One reason is that we have a branch 
point at threshold and the commutation of the derivative and the integral made above is 
not justified. Indeed, one can do the integration with respect to q, in Eq. (|10D, analytically 
and, after renormalization of the ultraviolet behaviour, the integral can be obtained from 



the two particle loop integral, Jo, given in appendix A of [32]. The difference is that in 
Eq. (pT0|) we use a non-relativistic propagator for the nucleon, but since we are concerned 
about the infrared divergence it does not make any difference. Then one can study the 
derivative with respect to k°. The expression for the loop function is: 

d 4 q 1 1 



t res {k) ~2MJ {k 2 ) = 2MiJ 



(2tt) 4 (q 2 -M 2 + ir]) ((q + k) 2 - m 2 + fry)' 
where the approximation sign is present because the l.h.s. of the equation is the baryon 
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non-relativistic approach to the r.h.s. expression. Removing constant subtraction terms 
irrelevant for the discussion, this expression |}2] is given by 



tA , 2t/ , m (M-m) 2 (7 + 1 M 2 — m 2 M 
(4tt) 2 J q (s) = (1 - -) a In + In 



(j—1 s m 



\ 



(M+m) 2 
s 



1 



(M-m) 2 ' 



(12) 



where s is the Mandelstam variable. Close to threshold (in this case s = (M + m) 2 ), and 
above it, we can expand Jq{s) as 

J o s =7r« S(-i7ra + 2a a + 0(^ + In— , 13 

(Any \ s s m J 

and then we can perform the derivative with respect to s and find (j^ = ^fz^yrp) at 
threshold from above: 

" w = - — -7- (-1 h2 In hC(cr) . (14) 

ds (4tt) 2 (M + m) 2 V 2a M + m m V V V ; 

Thus, what we find is that, since a goes to zero at threshold, the derivative of the imaginary 
part goes to infinity but the derivative of the real part remains finite. This is related to the 
fact that the imaginary part goes as the momentum of the particle, and then its derivative 
with respect to s is infinite at threshold. Below threshold a becomes purely imaginary 
and therefore dJo{s)/ds is purely real and it diverges at threshold. This is logical since 
the analytical continuation below threshold of the imaginary part above threshold becomes 
real and this is the reason for the infinite derivative. Hence, the derivative of Jo(s) with 
respect to s, and therefore that with respect to k° when k — 0, takes different values at 
the right hand side than at the left hand side of the threshold point. 

The results are then different from what one obtains by commuting the integral and 
the derivative, as we did before in Eq. fllPP to obtain the divergence of Eq. (pTTj) . However, 
the fact remains that there are divergences. Nevertheless, this analytical study has served 
to see that the origin of the divergence is the existence of the imaginary part in the free 
t-matrix close and above threshold. This realization is important because, as we shall 
see, for finite values of p, the effect of Pauli blocking drastically reduces the imaginary 
part of the rescattering term of the meson selfenergy close to threshold, and what is more 
important it can be differentiated with respect to k°. In order to see this, we perform first 
the integration over the Fermi sea in Eq. (|9|) and write the meson selfenergy as 

U-(k°,k =0,p) = J ^2^g) U (k°-u( q ), q ,p), (15) 

where U is the Lindhard function for forward going particle-hole excitation. The ordinary 
Lindhard function contains also the backward going particle-hole term. For this reason 
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we give in the appendix the explicit expressions for the forward going contribution to the 
Lindhard function needed here. Inspection of Eq. @ tells us that there is only imaginary 
part for k° — u(q) > 0. 

It is easy to prove, using the formula for ImU of the appendix, that, for k° — m << 
kp/2M, lmU res (k , 0, p) goes as (k° — m) 2 , but we save the proof here since this is a well 
known result for the imaginary part of the selfenergy of any particle in a medium |34| [36 . 
This result is important because then the imaginary part of the meson selfenergy and 
its derivative with respect to k° are continuous at threshold (actually for this rescattering 
contribution both the imaginary part of the function and its derivative vanish at threshold). 
Hence, the singularity in the derivative of the imaginary part of the free t matrix above 
threshold, which we have seen in Eq. ([14]), disappears in the meson selfenergy at finite 
densities. The use of dispersion relations in the selfenergy and its derivative then guarantee 
that both the real part of the selfenergy and its derivative are also continuous at threshold. 
Let us show this in a quantitative manner. 

In Fig. H] we show the real part of the meson selfenergy from this rescattering piece 
for k F = 100 MeV, corresponding to p = 0.05p , 50 MeV (p = 0.006p ) and 10 MeV 
(p = 0.00005po) divided by p for normalization, and we compare the results with the 
real part of the free t matrix for the meson nucleon interaction from the corresponding 
rescattering term, t res . We observe, indeed, that for finite densities the meson selfenergy is 
a regular function and shows no cusp, which is clearly visible in the free t res matrix. It is 
interesting to see that for very small densities (see the curve corresponding to kp = 10 MeV 
in the figure) U res /p goes like the free t res matrix, except for the fact that Pauli blocking 
has provided a regular ization around the cusp. It is also interesting to observe that even 
at kp = 50 MeV, and certainly at kp = 100 MeV, one obtains practically a linear function 
in k°. We can thus anticipate that problems linked to the original discontinuity of the 
derivative in the t matrix will have negligible effects in the evaluation of any observable 
where the different densities will have to be weighed in the nucleus. However, the formal 
problem still remains, because one has to define the derivatives at threshold for any density 
and, even if now they can be calculated, for very small densities (see e.g. the case of kp = 10 
MeV) the function around threshold is by no means linear in k°. Hence, if one takes a 
linear extrapolation based on the derivative at threshold, see Eq. (0), it is obvious from 
the figure that it would lead to a strong diversion from the actual calculated values of the 
selfenergy. 

As we have mentioned, the problem is a formal one affecting only to small densities 
which should have no relevance in practical problems provided a sensible and realistic 
approach is taken. Yet, this might not be always easy to achieve. As an example we 
could consider a sensible approach the one in fl30"| , where, according to the authors, they 
"artificially modify the slope of the function such that it is zero at /cf = but agrees to good 
accuracy with b e ff(kp) for kp > 50 MeV ". However, the final effect of the nonlocalities 
on the calculated level shifts and widths in that work is substantially larger than we find 
here when treating carefully the nucleus as a finite system, as described below. 

Indeed, the problems mentioned above do not appear in a finite nucleus. This issue was 
already addressed in J57] in a different problem, the study of the effect of core polarization 
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Figure 2: Real part of the kaon rescattering selfenergy, U res , of Eq. 
function of k° — m for different densities. 



divided by p as a 
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on inelastic or charge exchange pion nucleus reactions. The reason is simple but subtle. 
If one looks at the expression of the Lindhard function, which results from removing the 
integral over d 3 q/(2n) 3 , the factor l/2uj(q) and setting q = in Eq. and one takes the 
threshold value k° = m and k — 0, one realizes that there is an indetermination of the 
type 0/0 in the integrand. The analytical expression of the Lindhard function is such that 
its limit for k — > is finite. Yet, in a closed shell nucleus the numerator of the equivalent 
response function is zero because it involves a matrix element < <fri\exp(ikr )|0 2 > ; where 
<pi would correspond to an occupied state and (f) 2 to an excited state and these functions 
are orthogonal. However, the denominator would involve e« — ej, which is strictly non 
zero because there is a minimum excitation energy from the occupied states to the excited 
states. Thus, the corresponding response function in a finite nucleus is strictly zero while 
the Lindhard function does not vanish because of the continuity of the energies in the Fermi 
sea in infinite nuclear matter. Rather than doing an unnecessarily complicated evaluation 
in a finite nucleus, a simple solution to the problem was given in |J7|, including a finite 
excitation energy in the particle state, which we shall call the gap, A, and reevaluating 
the Lindhard function, which then turned out to be strictly zero in the limit of small k. 
The expressions for the Lindhard function with the gap can be found in the appendix 
of |37|, however they include forward and backward propagating ph excitations. We need 
here only the forward propagating ph excitation according to Fig. [I], and since we perform 
evaluations with this new function in what follows, we give the explicit expression for this 
function in the appendix. Once the gap is included, the rescattering term of the meson 
selfenergy, Il res , has imaginary part only for k° > m + A and, thus, around k° = m, the 
relevant region for kaonic atoms, it is purely real. The gap is physical and it is about 1-3 
MeV for most of the nuclei we analyze here. So, our strategy will be, first to show that once 
the gap is included the selfenergy can be accurately approximated by a linear function of 
k° around k° = m. Second, to test that the results for kaonic atoms are rather insensitive 
to the gap value, with changes on the shifts and widths far smaller than the experimental 
uncertainties. 

In Fig. § we show the results for the real part of the rescattering term of the meson 
selfenergy divided by p for different densities and compare them with the real part of the free 
t res matrix including the gap in the particle energy, this is to say, we increase the nucleon 
mass in the intermediate state of the rescattering term by the gap energy, i.e, we add (—A) 
in the denominator of Eq. (|T0|) . We have chosen here a gap of 3 MeV. The function around 
threshold is now well behaved for any density, showing a linear dependence in A; , and we 
can see that for values of kp around 50 MeV the meson selfenergy divided by p and the 
free t res matrix differ in less than two percent and they are practically indistinguishable at 
kp = 10 MeV, a consequence of the low density theorem for the meson selfenergy [38|-|0 . 

We did not yet pay attention to the dependence of the selfenergy on A; 2 , but the same 
arguments as above can be repeated. Again we can see in Fig. j| that, once the gap energy 
is considered, the dependence is smooth around threshold, k 2 = 0, and the selfenergy 
divided by p merges to the free t res matrix as a necessary consequence of the low density 
theorem. Here t res (k° = m, k) is defined as in Eq. (|10|) by changing q 2 /2M by (k — q) 2 /2M 
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Figure 4: Real part of the rescattering term contribution to the meson selfenergy, W es , 
divided by p as a function of k 2 for different densities, evaluated including a gap of 3 MeV 
in the nucleon particle energy. 



(and adding (—A) in the denominator). 
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The study conducted above has served to show the problems in the low density limit 
close to threshold due to the cusp behaviour of the free t matrix, but we have also seen 
that, once a physical value of the gap is introduced, all the problems disappear and one 
can also rely on the low density limit in order to evaluate values of the selfenergy and their 
derivatives with respect to k° and k 2 at low densities. 



4.2 Behaviour of the full kaon selfenergy 



Now we come back to the realistic situation of the kaon selfenergy including all coupled 
channels in the calculation. We have learned from the previous section that the kaon 
selfenergy does not behave linearly in k°, k 2 around thresholds. We have avoided these 
problems around the K~N threshold by means of the gap A. The rest of the thresholds, 
like 7rS, t]A, etc., are far from the region of energy and momenta met by the K~ atoms, and 
hence these other coupled channels do not affect the linear behaviour of the kaon selfenergy 
for kaonic atoms. We will obtain the low density limit simply taking the selfenergy equal 
to tp, where t stands now for the average of tx-N over protons and neutrons evaluated 
including the gap energy for the intermediate nucleon state in the Bethe-Salpeter equation. 
Thus, in the limit of p going to zero we obtain 



dU(k,p) 



dk 2 
dU(k,p) 



did 



(ui=mK, k=0) 



m 

ds 



M 
ds 



P 

(o) 

2(m K + M N )p 



(16) 



(0) 



where the label (0) indicates that the derivative is taken at threshold. This means that 
the values of b(p) and c(p) from Eqs. (H) can be cast in the low density limit as 



Kp) 



p 

c(p) 



P 



p=0 



p=0 



dt 

ds 



m 
ds 



(0) 

2{m K + M N ) 



(17) 



(o) 



and the above relationships are the constraints we impose to the low density behaviour of 
the parameters b(p) and c(p). Note that these coefficients are now complex because of the 
inclusion of the non-elastic channels. 

In Figs. [5] and ^ we can see the numerical results obtained for the partial derivatives of 
the kaon selfenergy as a function of the density. In order to stress the low density behavior 
we have also shown them divided by the density plotted against the Fermi momentum. 
The density and Fermi momentum are in units of po — 0.17 fm -3 and kp = 268.4 MeV, 
respectively. We also show the results calculated with the two values of the gap energy, 1 
MeV and 3 MeV. 

As one can see in the figures, the values of the partial derivatives have a density depen- 
dence that makes the derivatives positive and negative. Thus, there are subtle cancellations 
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when evaluating the results with the potential of Eq. (0) because, even if the derivatives 
act on the potential which depends on the nuclear density, the physical meaning is that one 
should be producing the kaon momentum which is a small quantity. Hence, the necessary 
cancellations have to appear in the integrals involving the derivatives of the density, which 
measure nucleon momenta, to finally account for the small kaon momenta. This is the 
reason why so much emphasis has to be made on the low density limit because, as one can 
see in Fig. [5], the low density limit for the partial derivative with respect to k 2 forces a 
change of sign of the derivatives with respect to p which will appear when the V operator 
acts on the p dependent functions of the potential. 



5 P— waves 



The lowest order meson-baryon Lagrangian |42| can be written as 



L [ B) = (Bi-fVpB)-M B (BB) 

+ ~Z>(B7"7s K> B}) + l -F(Bn,[u^ B}) , (18) 

where the symbol () stands for the trace of the SU(3) matrices B, involving the baryon 
fields of the octet of the nucleon, and u^, involving the meson fields of the octet of the 
pion. The constants D and F are given by D + F = qa = 1.257 and D — F = 0.33. 
The term involving the covariant derivative gives rise to a contact s-wave term, which is 
evaluated in |[24|| , and a contact p-wave term which we write here in Eq. (|19|) in the center 
of mass (CM) frame of the KN state 

t ~ Cl Up{ 2M % ) \ 2M j ) \2M { 2Mj J (19j 

with Mi j, Eij, the mass and energy of the initial and final baryons, and C^- SU(3) coeffi- 
cients which are tabulated in . The variables k 1 and k are the center of mass momenta 



of the outgoing and incoming meson, respectively. 

In addition we also have the contribution of the A, E and E* pole diagrams with two 
vertices of the type KNY, with Y being the hyperon, which come from the D and F terms 
of Eq. (pit) and that can be easily evaluated 

Unlike the A (1405) resonance, which is generated dynamically from the lowest order 
s-wave chiral Lagrangians and multiple scattering, the strength of the p-wave interaction 
is too small to generate dynamically the E*(1385), which, in the language of refs. |45|.[16|. 
would then qualify as a genuine, or preexisting, resonance built up mostly from three-quark 
states. Thus, we follow a phenomeno logical approach to include the contribution from E*/i 
excitations. The KNY,* vertex is evaluated in |PJ by means of SU(6) symmetry, in analogy 
to the evaluation of the 7rAA^ vertex from the ttNN one. It differs slightly from the one 



used in jKJ, where SU(3) arguments are used. 
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Figure 5: The real and imaginary parts of b versus pj po are depicted in the two left figures. 
The real and imaginary parts of [b/(p/po)} versus kpjkp^ = (p/po) 1 ^ 3 are shown on the 
right panels. The crosses correspond to the points which have been numerically evaluated 
from Eqs. (|1|), (|3p and ([[7]) and the lines are the interpolated values, which have been used 
for calculations. The solid lines have been obtained using an energy gap of A = 1 MeV 
and the dashed ones with A = 3 MeV. 
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Figure 6: The left figures show the real and imaginary parts of c versus p/po- The real 
and imaginary parts of [c/(p/po)) versus k F /k Fo = (p/Po) 1 ^ 3 are shown in the figures 
on the right-hand side. The crosses correspond to the points which have been numerically 
evaluated from Eqs. ([[]) , (§) and ([17}) and the lines are the interpolated values for arbitrary 
density c has been obtained from II calculated without pion selfenergy for pj p < 0.25. The 
solid lines have been obtained using an energy gap of A = 1 MeV and the dashed ones 
with A = 3 MeV. 
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For the case of A and £ pole terms the expressions that we get for the p-wave amplitude 
in the CM are 

where the subindices M, B (M', B') stand for the initial (final) meson, baryon. The quan- 
tities D MBY are SU(3) coefficients given by 

/20 D , — F , . 

Dm by = c DY \j y - c FY \/Yl— (21) 

where the cr> Y , cfy coefficients are given in table 0, and / = 1.15 f T as in ref. |24] . 
Similarly, in the case of the E* pole term the amplitude in the CM is given by 

= £We*£W- r \ f Sk' S^k , (22) 
Vs - Ms* 



where the D coefficient is given by 



_ 12 D + F 

U MBY,* — C S~^ — \ A6 J 



and cs is tabulated in table |l|. 



Table 1: Coefficients for the K-NY(K~N^*) couplings 



cda cfk cdt, c_fs 



cs 



K~p 



K~n 



W io 



In order to evaluate the p-wave selfenergy we write the amplitudes in terms of the kaon 
and nucleon variables in the frame where the nuclear Fermi sea is at rest. The former 
(1 + k°/2M B )ok vertex of Eq. ©, which was given in the CM, becomes now 



ok 1 



k 







- —k 

2M Y J 2 \M N ' My 



1 1 

+ 



(24) 



where p is the nucleon momentum, and the combination of the two vertices gives 



ok ok I 1 — 



2M 



+ 



Y 



1 



1 



M 



N 



M 



Y 



(25) 
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where we have already specified the forward direction of the kaons where the evaluation of 
the selfenergy in infinite matter is done. 

On the other hand, for the S* pole term we have to write the kaon momentum in the 
CM frame and then we have 

_ ( bO \ £0 

^H 1- (26) 

Thus, the kcukcu combination becomes 

kcM-kcM^kk{l-—j + (^— J f. (27) 

The evaluation of the selfenergy corresponding to the diagrams of Fig. [7] can be written 
for symmetrical nuclear matter in terms of the Lindhard functions used in |H]] as a p-wave 
part 

n«(*°,£,p) = ~D^ pA f A Pu A (k°Xp) 

+ \D 2 K ^ flk 2 U^Xp) 
+ D 2 K ^ fl,k 2 U^(k°,k,p) 

- lww N &2p (28) 

where the recoil factors are given by 



v 2M A J ' J \ 2M S ; 

h - = ( l -^)' h= ( l -M^) (29) 

together with an induced piece of s-wave nature coming from the Fermi motion of the 
nucleons and the frame transformation (p 2 terms in Eqs. fl25|) , (^7|) ) which is given by 

U^ nd \k°, k, p) = h 2 F (k°) 2 (A pPp + A nPn ) (30) 



where An(N = p, n) is given by 



4 = 7 7^ + 



1/1 1 \ 2 D 2 



K-NA 



4 \M N M A J y/S- M J 



+ 



4 \M N M^J y/s- M s 

2 

K-NH* 



( 1 \ 2 D 2 



~ ° K - N ' K - N IpM^ {m^Tw) (31) 
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Figure 7: Diagram contributing to p-wave K -selfenergy via K N — > A, K N — > E and 
K~N — > E* processes. 

and Ck- p ,k- p = 2, CK- n ,K-n = 1- For the kaonic atom case the Lindhard function U Y is 
given by: 



U Y (m K ,0,p) = — ,f -, y = A, E, E* (32) 

+ Mat - My + ^Fyfys = M N + m K ) 

Since the momentum in the k 2 factor of Eq. (|28| ) comes from the p-wave amplitudes 
and is already in the lab frame, the appropriate combination in finite nuclei, according to 
the findings of |T7|,]2|], is k 2 U Y -> - VC/ y V and k 2 p -> - VpV. The familiar ATT 
term (angular transform term) is here incorporated by means of the recoil factors f\ given 
in Eq. @. 

With all these ingredients we solve the KGE, Eq. (|6|), and show the results in the next 
section. 



6 Numerical results 

We solve the KGE of Eq. (||) with the microscopic antikaon-nucleus optical potential of 
ref. |13| (neglecting any isovector effects) plus the several non-local and p-wave terms 
discussed in the previous sections. As we already mentioned, for Vc(r) we use the Coulomb 
interaction taking exactly the finite-size distribution of the nucleus and adding to it the 



vacuum-polarization corrections |47 | 



The numerical solution of the Klein-Gordon equation is done using the method in 
coordinate space of fj8fl . The densities used throughout this work are those compiled in 
ref. [[J . However charge (neutron matter) densities do not correspond to proton (neutron) 
ones because of the finite size of the proton (neutron). We take that into account following 
the lines of ref. 
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In the present work we have used a K~ selfenergy calculated in symmetric nuclear 
matter. For heavy nuclei having a neutron excess our results could be easily improved 
by weighting the proton and neutron contributions to the selfenergy in Eq. (|l|) with the 
factors Z/(A/2) and N/(A/2), respectively. However, for the purpose of the present work, 
which is to establish the role of p-wave pieces and non-localities associated to s-wave, we 
have not considered necessary to implement these changes. 



6.1 P— Wave and non-local effects. 

We solve the KGE using different options to account for the non-local terms of the K~ 
selfenergy, EL For each one of these options, the KGE has been solved for a set of 63 
shifts and widths of kaonic atoms levels^. Shifts and widths for all models are given in 
table H (gap A = 1 MeV) and table |3] (gap A = 3 MeV) for a set of selected typical 
levels of kaonic atoms across the periodic table. To better quantify the changes we also 
give x 2 P er number of data. As can be appreciated in tables the effects on the spectra 
of both the p-wave pieces, discussed in section 5, and the non-local terms, coming from 
the dominant s-wave potential of ref. [13[ when the zero range approximation is relaxed, 
are small when compared to the experimental uncertainties. By comparing \ 2 /N = 3.76 
for row (1), where only the s-wave local optical potential is considered, with the values of 
X 2 /N for rows (2), (3a), (3b), (3c), (3d), (3e) and (4), where different nonlocal and p-wave 
effects are added, we see that x 2 /N change at most by ±0.7. Then we conclude that the 
effects of these nonlocal contributions on the kaonic atoms spectra are smaller than the 
experimental errors. Comparison of the results in the two tables shows that they are also 
rather independent of the value of the gap. Thus, despite some theoretical ambiguities in 
the nonlocal terms, for practical purposes they can be safely ignored. 

The small effect of the p-waves on the shifts and widths found here was also noted in 
ref. [|J, where the p-wave piece coming from E*/i excitations was included. Their p-wave 
contribution to the antikaon selfenergy turned out to be very small compared to the s-wave 
one (see Fig. 7 in ref. ||) and made only a "minor contribution" to the atomic shifts and 
widths. 



6.2 E*— hole excitation effects on the selfconsistent determina- 
tion of the s— wave K~— selfenergy. 

In this subsection we report on the role played by the E*/i— excitation when it is included 
in the self-consistent calculation of the s-wave K~ potential. In ref. only the Ah— 
and £/t— p-wave terms where included in the selfconsistent calculation of the dominant 
s-wave antikaon-nucleus potential. The results obtained with this potential were reported 
in refs. |§ and |§, and have served us here, see row (1) in tables |] and ||, as a reference 
to evaluate the effect of the non-local and p-wave terms studied in this paper. In a latter 



work, ref. | 44fl , in addition to the Ah— and £/i— p-wave terms, the £*-hole excitation term 



1 This set of data is the same used in 
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of Fig. [7] was included in the selfconsistent evaluation of s-wave K~— self energy, as well. 
With this new potential, a better description of the kaonic atom data is achieved. Indeed, 
this improved s-wave potential provides^ a x 2 /N = 2.89 to be compared to the value of 
3.76 provided by the original potential of ref. ||13|| . In Fig. |8| we show the s-wave K~ 
optical potential both, with and without YTh— excitation in the selfconsistent evaluation. 
Differences between both potentials are moderately small for the low densities relevant in 
kaonic atoms, but have a significant density dependence. At low densities, the inclusion of 
the S*/i— term leads to smaller, in absolute value, values for both the imaginary and real 
parts of the potential. Thus, there are two competing effects: reduction of the attraction 
because of a smaller real part and a reduction of the repulsion (increase of the attraction) 
because of a smaller imaginary part. It seems, that the latter effect is bigger than the 
former one and the resulting effect is an increase of the attraction which leads to a better 
description of the data in agreement with the findings of ref. 0. Thus, the new potential 
provides bigger widths and smaller, in absolute value, shifts than the one of ref. ||T3| . 

In the row labeled (l)s* of table ||, we present the results for selected kaonic atom levels 
which are obtained by using only the s-wave optical potential which includes the S* — h 
excitation in the selfconsistent evaluation of the s-wave K~ selfenergy, with a vanishing 
value of the gap. By comparing this row with row (1) of the same table, we see that 
the inclusion of the S* — h excitation in the selfconsistent evaluation of the s-wave K~ 
selfenergy improves the global agreement between the theoretically predicted values and 
the the empirical ones. We should note that unlike the nonlocal effects which depend 
on the prescription, the inclusion of the S* — h correction in the s-wave K~ selfenergy, 
although it has a moderate effect, is well defined. On the other hand, considering again 
only the s-wave K~ selfenergy, the effect of including non vanishing gap values of about 
A = 1 ~ 3 MeV has a very small effect on the shifts and widths of the known kaonic 
levels. This small effect can be quantified by observing that, for the considered set of data, 
the values y 2 /iV = 2.89 for A = 0, y 2 /iV = 2.83 for A = 1 MeV and y 2 /iV = 2.94 for 
A = 3 MeV, which we obtain using only the s-wave K~ selfenergy, are very close. 



7 Conclusion 

We have concentrated on the evaluation of corrections to the K~ nucleus optical potential 
originating from the momentum and energy dependence of the s-wave selfenergy, previ- 
ously developed, plus corrections to this potential originating from the p-wave part of 
the elementary K~N amplitudes. Some of the corrections lead to nonlocal terms in finite 
nuclei, and the fact that the calculations of the selfenergy are done previously using local 
approximations, like the assumption of a local Fermi sea, introduces some ambiguities in 
the choice of the form of the nonlocal potential. The use of several alternatives serves us 
to quantify the amount of intrinsic uncertainty in our approach. 

2 This value is obtained with a zero gap energy, and only very small changes are found for non- vanishing 
values of the gap energy. 
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X 2 /N 


10 

5 


B 


?3 7 M 


63 
29 

-64/ 


Cu 

r 4/ 


112 
48 


Cd 

r 5g 


238 
92 

-£7i 


U 

T 7i 


(1) 


3.76 


217 


551 


109 


368 


384 


1121 


528 


1437 


330 


1090 


(2) 


4.00 


213 


542 


110 


362 


392 


1110 


543 


1420 


350 


1076 


(3a) 


3.20 


211 


565 


102 


397 


361 


1229 


494 


1588 


302 


1291 


(3b) 


4.00 


234 


564 


118 


383 


415 


1172 


568 


1515 


357 


1196 


(3c) 


4.01 


234 


564 


118 


383 


415 


1173 


569 


1515 


358 


1197 


(3d) 


4.02 


233 


562 


118 


382 


414 


1170 


568 


1511 


356 


1194 


(3e) 


3.38 


219 


568 


110 


383 


391 


1182 


538 


1528 


336 


1203 


(4) 


3.69 


217 


552 


110 


371 


388 


1141 


534 


1465 


337 


1166 




2.89 


208 


575 


105 


398 


373 


1219 


512 


1550 


320 


1201 


Exp 




208 
±35 


810 
±100 


80 
±13 


443 
±22 


370 
±47 


1370 
±170 


400 
±100 


2010 
±440 


260 
±400 


1500 
±750 



Table 2: Widths and shifts of representative kaonic atom levels in eV obtained from different 
potentials, taking always the energy gap equal to A = 1 MeV. Different cases correspond to 
the local potential of ref. |]l3| (for a comparison in row (1) we give the results obtained with 
H(rriK, 0, p), first term on the r.h.s. of Eq. (|2|)) . Rows (2) to (4) correspond to different additions 
to this dominant piece: 

(2) Only p-wave Lindhard function non-local effects from Eq. (^g|) and the local induced term 
of Eq. @ are added. 

(3) Only bq 2 non-local effects, see Eq. are included using different ways: (3a) — V6V, (3b) 
-V6V - 0.5(A6), (3c) -6V 2 , (3d) b [(w - V c ? - p? (3e) 6Re[(w - V c ) 2 -fJ, 2 -U]. 

(4) Only c(uj — fi) "non- local" effects are added, see Eq. @. 

The results of row (l)s* ar e obtained from a purely s-wave optical potential, like in row (1), but 
including S* — h excitations as described in subsection 3.2. For each potential, x 2 P er number of 
data N is shown. We use a set of 63 shifts and widths previously used in ref. 1|. 
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351 
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(3e) 


3.25 
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386 
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531 


1548 


325 
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(4) 


3.65 


217 


552 


109 


371 


386 


1142 


531 


1479 


330 


1169 


Exp 




208 
±35 


810 
±100 


80 
±13 


443 
±22 


370 
±47 


1370 
±170 


400 
±100 


2010 
±440 


260 
±400 
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±750 



Table 3: Same as in table |2| but for an energy gap of A = 3 MeV 




<£ 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 ~ 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



p/po p/po 



Figure 8: The real and imaginary parts of the K~ optical potential, V op t = H/(2rriK), 
at energy k° = m K ~ and momentum k = versus p/po, are depicted in the figures. The 
selfenergy II is evaluated from Eq. (|l|). 

The solid lines have been obtained using no energy gap (A = 0), the dashed ones with 
A = 1 MeV and the dotted ones with A = 3 MeV. The dot-dashed lines correspond 
to the A = and omitting the S*-hole propagation of the K~ in the Bethe-Salpeter 
equations. 
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Fortunately, the corrections turn out to be small, smaller than the experimental errors, 
which gives support to the results obtained with only an s-wave optical potential, which 
has been the rule in the studies of kaonic atoms. We also make some choices of equivalent, 
energy dependent, potentials and find also small corrections. We have calculated widths 
and shifts for a large set of kaonic atoms. Although fits of better quality than the present 
results can be obtained, the potential has the merit of being a free parameter theoretical 
one which gives a fair reproduction of the data. 

We have made a thorough study of the optical potential around threshold and have 
shown the mathematical problems that one faces at very low densities. They are related to 
the cusp in the elementary K~N — > K~N scattering amplitude at threshold, where some 
of the derivatives become infinite. The peculiar low density behaviour found in |3(| is tied 
to these problems and a natural physical solution for them is found here by considering 
the finite excitation energy of real nuclei. 

As with respect to the use of empirical potentials to analyze the kaonic atom data, our 
results endorse the approaches which are based exclusively on an s-wave optical potential, 
once we see that all nonlocal terms generated lead to corrections which are smaller than 
present experimental errors. 

Another relevant finding of this paper is, that using a purely theoretical optical poten- 
tial, a quite satisfactory description of the kaonic atom data is achieved. For 63 data, we 
find x 2 /N = 2.89, using the s-wave local theoretical potential of ref. [P] which contains 
S* — h excitations in the selfconsistent calculation of the s-wave K selfenergy. 
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Appendix: the Lindhard function with a gap in the 
particle-hole excitation energy 

We define the Lindhard function for the forward going particle-hole excitation as 

U(q°, q, A, p) = 4 / nCk i P) [1 ^ + P)] . (A.l) 

where the energy gap A separates the occupied and unoccupied nucleon states. Through 
the paper in some situations we have referenced to U(q°,q,p) which is obtained from 
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Eq. ( |A.1| ) above, setting the gap to zero. We find for x < 2: 



ReU{q°,q,A,p) 



2Mk F 1 (x u-5 v-8 , 
1 1 In 



IX 2 2x [2 
4 \ x 



v — 5 + x 2 — 2x 



v-5 



x 



In 



v — 5 — x 2 — 2x 



v — 5 + x 2 — 2x 



(A.2) 



and for x > 2: 
ReU{q°,q,A,p) 

where 



2Mk F 1 I x u-5 1 
7T 2 2x I 2 2~^ + 2 



In 



z/ — <5 — x 2 — 2a; 



2Mq° 



q 2MA 

X ~V F ' "fcT' P 37T 2 



2k p 



For the imaginary part we find: 



Im£/(g°,g,A,p) 



Im Z7(g° - A, g, p) for g° > A 
for g° < A 



v — 5 — x 2 + 2x 
(A.3) 

(A.4) 



(A.5) 



where 



3 M 



lmU(q° jq , P ) = --Trp— (1 - z 2 )6(l- \ z |) - (1 - / 2 )#(1- | ^ |) 
4 qkp L 
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